EXCURSIONS OF DIFFUSION PROCESSES AND CONTINUED 

FRACTIONS 



ALAIN COMTET AND YVES TOURIGNY 



Abstract. It is well-known that the excursions of a one-dimensional diffusion 
process can be studied by considering a certain Riccati equation associated 
with the process. We show that, in many cases of interest, the Riccati equation 
can be solved in terms of an infinite continued fraction. We examine the 
probabilistic significance of the expansion. To illustrate our results, we discuss 
some examples of diffusions in deterministic and in random environments. 



1. Introduction 
Recently, Marklof et al. [3TJ |32] studied the random continued fraction 
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where the u„ are independent random variables with the same gamma distribution, 
i.e. for every Lebesgue-measurable set A C K+, 

(1.2) v{u,,eA) = j^^^^^y^-'e-y'^dy, /i > . 

An earlier work of Letac & Seshadri [22 had shown that, in the case of positive 
A, the continued fraction is a generalised inverse Gaussian random variable. One 
of the results obtained by Marklof et al. is an explicit formula for the probability 
density function for every complex A. 

The genesis of the work reported here was the curious observation that the same 
distribution — not the continued fraction itself — appears also in the articles of 
Bouchaud et al. [T^ (A < 0) and Kawazu & Tanaka [2J (A > 0) on diffusion 
in a Brownian environment with drift. More precisely, consider a process X with 
infinitesimal generator 
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where 



-dx 



and the functions a and b are the instantaneous variance and instantaneous drift 
respectively. Many fundamental quantities associated with the process can be ex- 
pressed in terms of the solutions of the equation 

(1.4) ^^H''/'(-,A) = A0(-,A). 

The order of the equation may be lowered by introducing the Riccati variable 

<P 

where the prime symbol denotes differentiation with respect to x. This yields the 
Riccati equation 

(1.5) ^ + + 2W'ix)U = 2X/a{x) . 
dx 

When the state space is M and W is Brownian motion with drift, this Riccati equa- 
tion (understood in the sense of Stratonovich [33]) admits a positive stationary 
solution whose law is the same as that of the continued fraction (jl.ip . This formu- 
lation in terms of the Riccati variable goes back to the pioneering work of Frisch & 
Lloyd |T9] and has since been used extensively, both in the physics [lOj |30] and the 
mathematics [371 HI] literatures. 

In this work, we show that, for many diffusions — be they in a fixed or a random 
environment W — a solution of this Riccati equation may be found in terms of the 
infinite continued fraction 

. . , ^ ^ 2X/a(x) 

(1.6) U{x,X) =Uo{x) + 



, 2X/aix) 

ui{x) 



U2ix) -\ 

Continued fractions arise in the study of many random processes on a discrete 
state space, notably birth-and-death processes (see for instance O [TU [181 EOl [23] 
and the references therein) and random walks in a random environment [H [J] 113) ; 
their occurence in the context of processes on a continuous state space is compar- 
atively rare, though not unknown [7]. In the remainder of this introduction, we 
provide a brief explanation of the origin and significance of the continued fraction 
in the context of diffusion processes, and summarise our main results. 



1.1. Probabilistic interpretation of the Riccati variable. The equation (II. 4|) . 
supplemented with the boundary conditions that may be required in order to specify 
the process X uniquely, has two non-negative solutions, say 0_(-,A) and 0-|-(-, A), 
that are of particular importance: they are characterised (up to a constant factor) 
by the fact that 0_(-, A) is non- decreasing and 0+(-, A) non-increasing. Let x denote 
the starting point of the diffusion and let y be any other point in the state space. 
Then the first hitting time H(jj) of y is the random variable defined by 



(1.7) 



H{y) :=inf{t>0 : Xt ^ y} 
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and we have the following well-known formula for its Laplace transform [9]: 

if X < y 



(1.8) (e'^^(^)) 



■ 0_(x.A) 

0-(a,A) 



if X > y 



We deduce 



This expresses two particular solutions of the Riccati equation — one positive (J7-) 
and one negative (J7+) — in simple probabilistic terms. The choice of sign in this 
notation is somewhat disconcerting, but its justification will soon become manifest. 

Following Pitman & Yor ^34j, one can gain further insight by defining a local 
time process L{x) at x by 

(1.10) Lt{x) -.^ lun ^ f 1(™+,)(X,) dr. 

This process is proportional to the time that X spends in the vicinity of the starting 
point up to time t. Let ( denote the lifetime of X. If X is recurrent (respectively 
transient), then L(^{x) is infinite (respectively finite) almost surely. The inverse 
local time L~^{x) at the starting point a; is the process defined by 

(1.11) L^^{x) -.^iniiT : Lr{x) > t} , 0<t<L^{x). 

It is a measure of the time that elapses before X has spent a total time t in the 
vicinity of its starting point. Using the trivial identity 

t= l(_oo,2:] (-'^r) dr + / l(^_oo) (-'^r) dr , 
JO Jo 

we can write 

(1.12) L~\x)=T_{x)+T+{x) 
where 

T-A^) = / ^i-oo.x] i^r) dr and T+,t(a;) = / l(a:,oo) i^r) dr. 



y=x± 



The processes TL(a;) and 7+ (a;) have obvious interpretations as occupation times; 
they are subordinators, i.e. non-decreasing Levy processes, and we have 

(1.13) E(exp[-Ar±,t(a;)]) = exp[-tV'±(a;,A)] , A>0, 

for some function ^p±{x, •) called the Laplace exponent of T±{x) [5 . By the Levy- 
Khintchine formula, the Laplace exponent can be written in the form 

(1.14) i;±{x,X)= f (l-e-^y) iy±{x,dy) 

J(0,oo] 

for some measure i'±{x) on (0, oo] such that 

/ min{l,y}i/±(a;,dy) < cx). 

J(0,oo) 

The measure i'±{x, •) is called the Levy measure of the process T± (x); it describes 
the distribution of the heights of its jumps. Each jump corresponds to an excursion 
of the process X, i.e. to the path followed by the process between two successive 
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visits to the starting point x, and the jump height is the duration of the excursion. 
The + (respectively — ) case picks out the "upward" (respectively "downward") 
excursions, i.e. those whose path hes entirely above (respectively below or at) x. In 
the transient case, one or both of v+{x^ •) and i/_(x, •) will have an atom at infinity, 
reflecting the fact that X begins an excursion of infinite duration at the time of its 
last visit to x. Pitman & Yor [34] show that 



(1.15) 



i^±{x,X) = :f^U±{x,X) 



Knowing U±, one can therefore recover the Levy measure i'±{x, •) by inverting a 
Laplace transform. 

1.2. Stieltjes functions. Knight ^2E^ and Kotani & Watanabe [SQ make the ob- 
servation that Krein's theory of strings implies the existence of a measure, say 
<7±{x), on [0,oo) such that 

(1.16) -±(^-dz) 



C/±(a;,A) -T2A / 



'[0,oo) X + Z 

So, knowing U± , we can obtain a± (and hence i>± ) from 



(1.17) 

and the Stieltjes-Perron inversion formula 



(T±(a;,{0}) = ^C/±(x,0) 



TU±{x, —A — ie) 



dA 



(1.18) a±ix,A)= lim — I Im 
which holds for every ct-i- -measurable set A C (0, oo). Then 

poo 

(1.19) {x, {oo}) = a± (x, {0}) and i'±{x,dy)= / ze^*'^(T± (x, dz) dy . 

cr(dz) 



A function of the form 

SiLu) = 



UJ € 



UJ + z 

is called a Stieltjes function (or transform) . Stieltjes [3 7) considered the problem of 
recovering a measure on from its moments, and the continued fraction 

(1.20) ^ 



1 



niQUj 



1 



miuj 



£2 + • • • 

played an important part in his solution of this moment problem, as we proceed to 
explain. Suppose that the measure cr satisfies the following moment condition: 



(1.21) 



POO 

For every non-negative integer n, / z"cr(dz) exists . 

Jo 



For simplicity, suppose also that a has infinitely many points of growth. Then S 
has an asymptotic expansion in decreasing powers of uj, namely 

r poo 



-l-oo . 
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and one may associate with S two sequences {m„} and {in} of positive numbers 
by requiring that the finite truncations of (|1.20p have asymptotic expansions that 
agree with the above series up to some order ^] [3]. Importantly, only the first 
2n + 1 of the moments are required in order to calculate the coefficients rn„ and 
in- Stieltjes showed how, knowing these coefficients, a measure with the required 
moments can be constructed [Tl|33l[37]. Furthermore, he showed that if the series 

oo CO 

rrin or ^ in 

n—1 n—0 

diverges, then there is only one measure cr with the given moments, and the infinite 
continued fraction converges to S{uj). 

1.3. Statement of the main results and outline of the paper. Our first task 
is to bring out the relationship between the continued fraction of Stieltjes and our 
own. To this end, and with a view to making the paper reasonably self-contained, we 
recall in fj2lsome of the results from Krein's theory of strings that are most relevant 
to diffusion processes. In Krein's terminology. Equation (|1.20|) is the characteristic 
function of a string with a discrete distribution of masses to„ such that mo is at 
and the spacing between m„ and m„_i is £„. Then, in ^ we prove the 

Theorem 1. Let X be a non-singular diffusion process started at x. Suppose 
that (T±{x, •) has infinitely many points of growth and that the following moment 
condition holds: 

(M) For every positive integer n, / z~"'a±{x,dz) < oo . 

Jo 

Then there is a sequence of positive numbers 

Tui,±{x), TU2,±{x), . . . 
such that, for every, n — 0, 1, 2, ... 

C/± (X, A) - C/± (X, 0) = , , + 0(A"+i) 

2A/a(x) 
ui,±[x) H 



, 2X/a{x) 
U2.±[x) H 1 -—- 

Un,±[X) 



as \^ 0+. 

Conversely, suppose that the identity 



2X/a(x 

U±{x,X)-U±{x,0) = ' ^ 



_ 2X/a{a 

Ul,±[X) 



U2,±ix) H 

holds for some infinite sequence Tui,±{x), T'"2,±(a;), . . . of positive numbers such 
that the series 

oo 

71=1 

diverges. Then a±{x, •) has infinitely many points of growth and the moment con- 
dition (M\) holds. 
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By using ()1.19p and Fubini's theorem, it is easy to see that the moment condition 
M|) can be written in the equivalent form 



(1.22) For every positive integer n, / y"^ i/±{x, dy) < oo . 

Jo 

So when the coefficients are positive, we immediately glean some information about 
the tail behaviour of the Levy measure. In this case, the continued fraction (jl.6p 
is related to that used by Stieltjes in his study of the moment problem, and the 
Riccati solution L'±(a;, A) has an expansion in increasing powers of A, from which 
the continued fraction coefficients can be computed by a well-known algorithm [3] . 

The moment condition ((M)) . or its equivalent form (|1.22p . is by no means neces- 
sary for the Riccati solution U± to have a continued fraction expansion if we allow 
the coefficients Un,± to be of arbitrary sign. In such cases, U±{x,X) is no longer 
asymptotic to a power series in A. Nevertheless, the continued fraction coefficients 
may be computed by an algorithm, presented in a more general form in Common 
& Roberts [H] , which we describe in 21 

Theorem 2. Suppose that the Riccati equation il.5]) has a solution expressible in 
the form m.6\) . Then, for n g N, u„ solves the homogeneous Riccati equation 

(1.23) ^+ul + 2W^ix)un = 0, 

wh- 



ere 



(1.24) w'^^^^-u^^,^W'^_,, W'^-.^W . 



In some cases, this expansion algorithm admits a probabilistic interpretation in 
terms of a sequence of diffusion processes, where the nth process corresponds to 
the environment Wn- In ^ we study the relationship between adjacent processes 
in this sequence. This reveals a surprising connection with the Ciesielski-Taylor 
theorem generalised by Biane [H [12] , which exhibits a large class of ordered pairs 
of diffusions such that some occupation time of the first process has the same law 
as some hitting time of the second process. 

In fJH we consider some simple examples of diffusions in a deterministic environ- 
ment which illustrate these results. 

Finally, in ^ we return to the problem that provided the initial motivation 
for this work, and consider the case where the environment W is itself a diffusion 
process: 

Theorem 3. Let W be a standard Brownian motion with positive drift. Then the 
system of Stratonovich equations \1.23'l[T724\ l has a stationary solution such that 
uq = and, for n > 1, the Un are independent random variables with the same 
gamma distribution. 

Thus, we recover the continued fraction of [251 ISIl 131] • 

As this outline indicates, the present paper deals exclusively with what may be 
called the direct problem: given the characteristics a and b of the diffusion, together 
with appropriate boundary conditions, we use a continued fraction expansion to 
compute the Levy measure of the excursions. The inverse problem, namely that 
of finding the diffusion given the Levy measure, is also of great interest [15]. The 
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relevance of our results to the solution of this inverse problem will be developed in 
a separate publication. 

2. KrEIN'S theory of strings and DIFFUSIONS 

This section provides a succinct review of some well-known results concerning 
strings and diffusions. Kotani and Watanabe's account of these topics in ^8] is 
particularly well-suited to our purpose, and we follow them very closely. The reader 
familiar with this material need only take note of Definition 12.21 before proceeding 
to the next section. 

2.1. Strings. A string, say m, is a function from [0, oo] to [0,oo] that is non- 
decreasing, right-continuous and infinite at infinity. Let m be a string and suppose 
that m is not identically infinite. Set 

c :— inf {x : m(x) > 0} and £ sup{x : m(x) < oo} . 

The number £ is called the length of the string. 

By setting m(0— ) = 0, we obtain from m a (Stieltjes) measure dm on [0,cx)). 
Let uj > and denote by and ri[-,uj) the (unique!) solutions of the following 

integral equations on [0,^): 

y+ 

^{z,uj) dm(z) 



^(x, cj) = 1 + 0; 
and 

77(x, Oj) — X + LJ 



dy 



y+ 



dy. 



^ (respectively rj) can be viewed as the particular solution of the generalised differ- 
ential equation 

d^ 

(2.1) __/(., ^)=^/(.,^), 0<x<£, 

subject to the condition 

/(0,w) = 1 and ^(0,w) = m(O-l-) w ^respectively f{0,uj) = and ^(0,^) = 1 

We shall refer to (^, 77) as the pair of fundamental solutions associated with the 
string. 

The function 

S(uj) := hm — 

x^£-e(x,a;) 

is called the characteristic function of the string. It may be shown that S is neces- 
sarily of the form 

cr(dz) 



(2.2) Siuj) = c + 



[0,00) 



UJ + z 



for some measure a on [0, 00) such that 

cr(dy) 



< 00 . 



'[0,00) 1 + ?/ 

Conversely, Krein showed that every function of the form ()2.2p is the characteristic 
function of a unique string. 
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The right-continuous inverse, denoted m* , of a string m is called the dual string 
of m. Denote by 

J[0,oo) 

the characteristic function of the dual string. Then 

c* = m(0+) and S*{uj) = J. . 

2.2. Generalised diffusion processes. Given a pair (m_,m_|_) of strings such 
that m_(0+) = and £± > 0, let din_ be the image measure of dm_ under the 
map X I— >■ — X and set 

dm_ on (— oo, 0) 
dm+ on [0,oo) 

The support of dm, denoted Im, is a subset of {—£-,£+). Let _B be a standard 
Brownian motion, denote by Lf{x.) its local time and set 



dm := 



ip{t) :-^Lf(x) dm. 



Then 

defines a Markov process on /m whose lifetime is the first hitting time of — ^_ or 
£+. X is called the generalised diffusion corresponding to the pair (m_,m_|_). 

Definition 2.1. A one-sided diffusion is a generalised diffusion such that either 
m_ or m+ is the zero string. 

Thus a generalised diffusion may be thought of as an ordered pair of one-sided 
diffusions. 

The following non-standard definition will also be helpful later on: 

Definition 2.2. The Krein dual X* of the generalised diffusion process X corre- 
sponding to the pair (m_,m+) is the generalised diffusion process corresponding 
to the pair (ml, m!j_). 

2.3. Non-singular diffusion processes. Next, consider a diffusion process X 
with generator (|1.3|) whose state space / is an interval with left endpoint I and 
right endpoint r. For the sake of greater clarity, we shall, in this subsection, depart 
from our usual notation and use Xq instead of x to denote the starting point of 
the process; this leaves us free to use x to denote some generic point in the closure 
of /. We say that X is non-singular [26l [28] or regular [9] if, for every xq, y G I, 
H{y) is finite with positive probability. We shall show that such a process is also a 
generalised diffusion process in the sense of the previous subsection. 

Associated with the diffusion is a scale function s and a speed measure m given 
respectively by 

s{x) = / e-^^'^y^ dy and m{dx) = -e^^ dx . 
Jxo a 

The infinitesimal generator may then by expressed in the equivalent form 

d / 1 d 



dm \s'{x) dx 



I < X < r . 
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In general, the scale function and the speed measure are not enough to determine 
completely (up to a constant factor) the non-negative monotonic solutions </>_ and 
0+ of Equation (|1.4p that characterise the process; some boundary conditions must 
also be imposed. Following the accepted terminology, we say that r (respectively 
I) is entrance if 



Xa 



s{x) ~ s{xo)] m(dx) ( respectively / [s{xo) — s{x)] m{dx) I < oo 



xo 



and that r (respectively /) is exit if 



m{{xo,x)) s' (x) dx respectively / m{{x,XQ))s' (x) dx \ < oo 



Xq 



A point is called non-singular if it is both entrance and exit, and natural if it is 
neither. A boundary condition for </)+(•, A) (respectively (/)_(•, A)) is required only 
if r (respectively I) is non-singular, and then only at r (respectively 

(1) If r (respectively I) ^ I, then this Feller- type condition is 

a+0+(r. A) + /3+s'(r) </)^(r. A) + 7+A<^+(r, A) = , 

(respectively «_</)_(;, A) /3_s'(0 0'_(/, A) 7_.A0„(/, A) =0,) 
where 

a± > , /3± > , 7± > and a± + /3± + 7± = 1 . 

(2) If r (respectively I) ^ I then the killing condition is 

(/)+(7-— ,A) (respectively (Z+, A)) =0. 

Next, we explain how to construct a pair of strings associated with X. Set 
X = s{x). Then 

= -i—r ' 1< X < r , 
dmdx 

where 

m(dx) = m(da;) , ^ ■— s'{y)dy and ^ s'{y)dy. 

J Xo Jxo 

Let m± be as follows: for < x < r (respectively < x < —1), set 

m+(x) =m([0,x]) (respectively m_ (x) =m([— x, 0])) . 

To make m± into a string, we need to define it over the whole of [0,cx)]. It will 
suffice here to discuss the extension of m+ into a string of length ; the extension 
of m_ into a string of length ^„ is analogous and will be immediately obvious. 
If r is singular, then it must be that r + m ([0, r]) = oo. This forces :— r and 

m+ (x) = oo for x > ^4. . 

On the other hand, if r is non-singular, we need to bring in the boundary condition 
satisfied by 0+(-, A): 

(1) If r belongs to /, we use the Feller-type boundary condition, i.e. we set 
:= r + ^ and 



m+(x) 



m([0,r])-H^ forr<x<^4 
00 for X > 
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(2) If r docs not belong to /, we use the killing boundary condition, i.e. we set 
= r and 

m+(x) = CXI if X > . 

Now, denote by , ri± ) the pair of fundamental solutions corresponding to the 
string m± thus defined, and by S± its characteristic function. Set 



e(x,A) 



^_(-x,A) if-^_<x<0 



and 77(x, A) 



-r;_(-x,A) if-^_<x<0 
77+(x,A) ifO<x<^+ 



^^+(x,A) ifO<x<^ 

Then the non-negative monotonic solution 4'±{-, A) of Equation (|1.4|) . suitably nor- 
malised, is obtained by setting 

1 



b±{x,\) = C(x, A) =F 



^±(A) 



77(x, A) 



and we deduce easily that 
(2.3) 



s±i\) 

It is clear from the foregoing construction that Cj_ = 0; so this formula is equivalent 
to our earlier equation (|1.16p . from which (|1.19l) follows immediately. 

Remark 2.1. In the case where the process X is started at a reflecting boundary, 
then one of m+ or m is the zero string, and so the corresponding generalised 
diffusion process is a one-sided diffusion. 



3. Proof of Theorem 1 

Proof. Suppose that cr±{x,-) has infinitely many points of growth and that the 
moment condition ((M| holds. Define a measure a on by 



(3.1) 



aiA) 



^cr±(a;,dz) , A -.^ {1/ z : z E A} 



Then a has infinitely many points of growth and satisfies the Stieltjes moment 
condition (|L2T|) . So we can write 



a{dz) 
UJ + z 



ajdz) 
1 -I- z/uj 



E 

3=0 



{-zya{dz) 



oj as w — > -l-oo . 



Following Stieltjes [33 (see also [U [551 133]), one can construct from this series two 
sequences {mfc} and {t-k} of positive numbers such that, for every n S N, 



E 

i=0 



(-z)V(dz) 



0(w-"'-2) 



as uj 



-(X), where 



2n - 1 if r„ = 

2n if r„ — m„uj 
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Thus, in the same limit, we can write 



1 

mooj H 



1 

h + 



1 

miLj ^ 



1 

£2 + ■ ■ ■ + — 



z-i.1/2 



Jo 1 + Jo 1/'^ + 

where we have made use of Equation (|1.16p to obtain the last equality. The first 
statement in the theorem follows if we take, for n = 0, 1, . . ., 

(3.2) ^ U2n+l,±{x) = and ^ U2n+2,±{x) = 2 in . 

a[x) 

To prove the second statement, construct from the given Wn,±(x) two sequences 
{to„} and {in} of positive numbers via Equation p.2p . The hypothesis implies 
that the expansion 

1 



1 



niQUj 



1 



1 

miLj 



i2 + --- 

is a well-defined function of cj G C\IR_. As shown by Stieltjes ^7\, since 



00 



^ TO„ or ^ in 

n=0 n=l 

diverges, there is one and only one measure a on such that 

a{dz) 1 

UJ + z 1 
mow H 



1 



i2 



Furthermore, a satisfies the moment condition (jl.2ip . The uniqueness of a implies 
easily that Equation p.ip holds, and so (J±{x, ■) satisfies the condition (jM]) . The 
fact that it has infinitely many points of growth follows from the fact that the 
sequence of the Un,±{x) is infinite. □ 



4. The expansion algorithm 

We now describe a method of obtaining the continued fraction which is based 
on the near-invariance of the Riccati equation under a certain linear fractional 
transformation. The same idea, which can be traced back to Euler [17], has been 
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used and developed more recently in a different context by Common & Roberts 
[m. Set 



(4.1) 



U{x, A) =: Uq{x, A) = uo{x) + 



2X/a{x) 
Ui{x,X) ' 



where uq and Ui are some functions which we shall specify presently. Substitution 
in Equation (|1.5I) yields 



2X/aa' 2X/a , 



2X/a 



Un 772"^i + "o + 2mo^^77^ + 2A/a 

^1 



2X/a 

2W'uo + 2W'^^^2X/a. 



Ui 



Now choose uq so that it solves the homogeneous Riccati equation 



ul + 2W'uo = . 



Then the equation satisfied by [/i is 

U[ + ^1 + 2 - "0 - W'^ Ul = 2X1 a . 

This is of the same form as Equation (|1.5p . save that W has been replaced by 
a! I (2a) — uq — W' . By iterating, we deduce the expansion 



(4.2) 



U{x, A) = uo{x) + 



2X/a{x) 



ui{x) 



2X/a{x) 



Unix) + 



2X/a{a 



Un+i{x, A) 

In this expression, m„ satisfies the homogeneous Riccati equation 



(4.3) 



^+ul + 2W:,ix)u„^0, 



the remainder Un satisfies the inhomogeneous Riccati equation 



(4.4) 
and 



da; 



+ Ui + 2W;,ix)Un = 2X/c 



Theorem [5] is thus proved. 

The homogeneous Riccati equation always admits the trivial solution. The cal- 
culation of the non-trivial solutions is sometimes facilitated by the following 

Proposition 4.1. The nontrivial solution u„ of the homogeneous Riccati equation 
of index n satisfies the recurrence formula: for every n > k > 0, 



Unix) = 3- In 
da; 



-2Wk(x) 
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Proof. Since, by assumption, uj ^ for fc < j < n, we deduce from the homoge- 
neous equation that 



dx 



Hence, for k < j < n, 

Uj{x) 

Let us write 



-2Wi 



/ e-2W'^ da; da; 



Wn{x) = / W:,(x)dx, 



-2^^ dx 



where the indefinite integral means that Wn is any primitive of the integrand. For 

n > k, 



2Wn = -2 y" W^ix) dx = -2 J 



2~a 



un-i - w;_i 



da; 



hia + 2M^„_i + 21n 



dx 



and so 



dx 



2 ^-2W^-^ 



The required formula follows easily by iterating. 



□ 



5. A PROBABILISTIC INTERPRETATION OF THE ALGORITHM 

It is straightforward to give a probabilistic interpretation of the first few coeffi- 
cients in the continued fraction. Indeed, from Equations (jl. 14111. l"5|) and Equation 
(15) of [31], 

uo^_(x) - uo,+ (x) 

is inversely proportional to the mean of the local time spent at x. Also, 

1 

=Fa(x)ui,±(x) 

and hence =pMi -|-(x) is inversely proportional to the average duration of the finite 
excursions to the right or left (— ). 

Our aim in this section is to gain some insight into the probabilistic content 
of the expansion algorithm itself. If the initial diffusion X has a Riccati variable 
U with a continued fraction expansion, then, as can be seen from the proof of 
Theorem [2] given in the previous section, the algorithm produces a sequence of Wn 
and a sequence of remainders [/„ . Each J7„ can be thought of as a Riccati variable 
of some diffusion, say A"„, corresponding to the "environment" Wm and it is then 
natural to investigate the relationship between adjacent diffusions in this sequence. 
In pursuing this line of thought, it is important to bear in mind the following points: 

(1) The Wn depend on which of the two Riccati variables Uj^ and U- is being 
expanded. In particular, if both U+ and U- have a continued fraction 
expansion, then there are two sequences of Wn ■ 



yv±{x,dy) 
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(2) It takes two strings to specify a regular diffusion — eacfi string corresponding 
to a one-sided diffusion (see Definition 12.11 and Remark l2.ip . Each Riccati 
variable specifies one string, and the regular diffusion is obtained by welding 
the strings together. 

In the remainder of this section, we will show how, from one continued fraction 
expansion of a Riccati variable of the regular diffusion X, one can construct a 
particular sequence of diffusions Xn such that the remainder Un is a Riccati variable 
of Xn- It should be clear from the points just made that sequences with this 
property cannot be unique. Our particular construction is inspired by the following 
key formula, proved by Pitman & Yor |34| : 

(5.1) ^±{x, A) = tp±{x, 0)T~ E„ (e-^«(-) Hix) < oo) . 

To simplify matters, we assume that the speed measure m of X is absolutely con- 
tinuous with respect to the Lebesgue measure and, with some abuse of notation, 
write 

m(A) — / m(x) dx 

J A 

for every measurable set A. The instantaneous variance a and the instantaneous 
drift b oi X are then given by 

a = and b — . 

ms m dx s 

5.1. The /i-transform. For definiteness, let us suppose that we work with a con- 
tinued fraction expansion of the Riccati variable [/+ . Taking the -|- sign in Equation 
(|5.ip and making use of the identity (|1.15p . we obtain 

U+{x, A) = U+{x, 0)-4- fe-^^(^) H{x) < oo) . 

dy y=x+ V / 

We will show that the second term on the right-hand side is the Riccati variable 
of a diffusion, say Y, obtained from X by conditioning, and we will express the 
characteristics of Y in terms of the characteristics of X. In particular, it will follow 
from our construction that the Riccati variables of Y are related to those of X 
via 

(5.2) U+{x,X) ^U+{x,0) + UX{x,X) and U-{x,X) ^U+{x,0) + Ul{x,X) . 

If [/+(•, 0) vanishes, then we take Y — X. Otherwise, X is transient and the 
event 

lim Xt = r , 

where C denotes the lifetime of X, occurs with positive probability. Let I < z < x 
and denote by Y the process obtained from X by conditioning on the event 

H{z) < oo . 

We have (see [9], 11.12) 

P:, {H{z) <oo)= (l)+{x, 0) h{x) , 

where P^; is the probability measure associated with X started at x and (f>+ is 
normalised so that (/)+(2:,0) = 1. 
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^ /i-transform ^ Krcin duality _ , , 

A > Y > ) 

s' h~'^s' K^m 

a a a 

b b + ah'/h a'/2~b - ah'/h 

Table 1. The transformation 3^ and its effect on the character- 
istics of the diffusion. 



For z < y < X, denote by dy an interval of infinitesimal length centered on y. 
Then 

P. (r. e dy) = p. (x, e dy h[z) < oo) = ^ < 



(i/(z) < (X) 
{Xt e dy) Py (i/(z) < oo) _ P, (Xt e dy) h{y) 



Px (iJ(z) < oo) h{x) 

Since ft- is excessive (see [9], 11.30), this calculation shows that Y is the ft-transform 
of X. This result does not depend on the particular choice of z as long asl < z < x. 
Clearly, 

The characteristics of Y are displayed in Table [TJ 

The foregoing discussion assumed that one is working with the Riccati variable 
U+. If, instead, one is working with the other Riccati variable J7_, then one should 
use h = ,\) to define Y . It is then easily verified that 

U+{x,X) = U^{x,Q) + Ul{x,X) and U^{x,X) = U^{x,Q) + Ul{x,X) 

holds instead of Equation 



5.2. Krein duality. Given Equation (|5.2p . we ask next for a diffusion, say Z, 
whose Riccati variables are related to those of Y via 

(5.3) C/|(x, A) Ul{x, X) = 2X/a{x) . 

We emphasise that this equation should be understood as specifying both Uf_ and 
Equation (|2.3p suggests that the Krein dual of Y might be a suitable can- 
didate. This raises some technical issues: F is a well-defined diffusion process, 
but its behaviour at the boundary points may differ from that of X, and it is not 
immediately clear whether Y is non-singular or whether it is even a generalised 
diffusion process in the sense of ' i'2.'2\ To avoid these complications, we shall be 
content to suppose that Y belongs to the particular class of non-singular diffusions 
described in the statement of the following proposition, so that its Krein dual exists 
and is a non-singular diffusion in the same class. The resulting concept of duality 
is then analogous to that used by Jansons [35], Soucaliuc [3S] and Toth [3B] (see 
the definition of conjugate diffusion in his Appendix 1). 

Proposition 5.1. Let X be a non-singular diffusion process, started at x, on an 
interval I that includes no non-singular boundary points — except possibly reflecting 
ones. Let X* be the Krein dual of X, in the same scale as X. Then 
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(1) X* is a non-singular dijfusion process, started at x, on an interval I* 
with the same endpoints as I and that includes no non-singular boundary 
points — except possibly reflecting ones; 

(2) the speed measure of X is the scale function of X* and vice-versa; 

(3) the Riccati variables and Ul_ of X* are related to the Riccati variables 
U+ and U- of X via 



where a is the common infinitesimal variance of X and X* . 

Proof. We work "backwards", i.e. from X, we define a non-singular process X on 
an interval /* with the same endpoints as /, started at x. Then we show that its 
associated pair of strings is the same as that of the generalised diffusion process 
X* . For the sake of convenience, we shall use in the proof the same notation for 
the process X as in ^ 12. 31 

In the interior of /* , the process X is determined by the generator 



The boundary behaviour of X is as follows: for every non-singular p G {1,t}, if p 
is reflecting for X then it is killing for X and p ^ /*; if p is killing for X, then p is 
reflecting for X and p £ I* . 

We remark that, for 1 < x < r, 



where dm(x) — m(dx) = m(dx). Recalling the definition of the string m±, this 
shows that the natural scale of X is 



and we deduce 

(5.4) m+(y) = m;(y) , < y < m+(r) , m_(y) - ml(y) , < y < m_(-l) . 

To show that X — X* , there only remains to extend the validity of these two 
equalities to every y e [0,cx)]. Consider the extension of the first equality. As 
explained in §2.31 the definition of m+ in [m+ (r) , oo] depends on the nature of the 
boundary point r. We bear in mind that, since the speed measure of X is the scale 
function of X and vice- versa, r is entrance for X if and only if it is exit for X and 
vice-versa. There are three cases to consider: 

(i) If r is singular for X, then m_|_(y) = oo for y > m+(r). Since r is also singular 
for X, we have r -I- m+(r) = oo. If m+{r) — oo, there is nothing to extend; if 
m+(r) < 00, then r = oo and, since is by definition the right-continuous 
inverse of m-|_. 



U±{x,X)U;^{x,X) =2X/a{x), 






m+(y) = oo = m+(y) for y > m+(r) 
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(ii) If r is reflecting for X, then 

r for m-|_(r) < y < oo 



m+(y) 



oo for y — oo 



On the other hand, r is killing for X, i.e. m_)-(x) = oo for x > r. Again, since 
m!j_ is by definition the right-continuous inverse of m+, we deduce m^(y) = 
m+(y) for y > m+(r). 
(iii) If r is killing for X, then it is reflecting for X, and the desired result follows 
from (ii) by symmetry. 

The extension of the second equality in (|5.4[) follows along the same lines. We have 
thus shown that X = X*. 

The first two statements in the proposition follow immediately. For the last 
statement, we apply Equation (|2.3p to X and to X*; this gives 

s'{x) m(a;) 2/a{x) 



U±ix,X)Ulix,X) 



S±{X)Sl{X) 1/A ■ 

□ 



Table [T] shows how the characteristics of a process transform under Krein duality 
in the situation envisaged by Proposition 15.11 

To summarise the foregoing discussion, one can — at least in some cases — de- 
scribe the first iteration of the expansion algorithm in terms of a map ^ obtained 
by composing two transformations: 

^ /i-transform , , Kroin duality a?- / v\ 

A >■ ¥ )■ Z =: ,yii (A j . 

More precisely, expanding the Riccati variable U+ of X leads, after one step, to the 
Riccati variable of the "new" diffusion ^h{X) where h — 0+(-,O). 

If, instead, one works with the Riccati variable U-, then h ~ </>-(•, 0) is used to 
construct the next diffusion 5^ (X). 

5.3. The Ciesielski— Taylor theorem. Ciesielski and Taylor [T^] noticed that the 
total time spent by a ((i + 2)-dimensional standard Brownian motion inside the unit 
ball in K'^^^ has the same distribution as the first hitting time of the unit ball in 
M'^ by a d-dimensional standard Brownian motion. This result can be expressed in 
terms of one-dimensional diffusions as an identity in law between some occupation 
time of a Bessel process (of parameter p) and a hitting time of another Bessel 
process (of parameter 1). Biane [6^ generalised the Ciesielski-Taylor theorem to 
other pairs of diffusions. As we shall see presently, his construction of these pairs 
can be described in terms of the map defined in the previous section. A similar 
observation has already been made by Toth in [38 . 

Suppose that (1) s{l) = — oo; (2) s(r) = 0; (3) if r < oo, r is killing, and 



s{y) m{y) > if Z > — c 



(4) 



y- 



,/^oo = °° ifZ = -oo 

The hypothesis ensures that 

Xt > r almost surely . 
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Then h :~ 0+(-,O) = — s > 0, and the /i-transform y of X is a non-smgular 
diffusion such that r is singular and, if I is non-singular, then it is killing. Its dual 
Y* = '%l{X) is therefore a well-defined diffusion. 

We seek a diffusion Z, with a suitable initial law, such that, for every I < x < 

y <r, 

K 



(law) 



(5.5) inf{i>0 : Zt ^ y} 

Following Biane [S] , the Feynman-Kac formula gives 



^{Xt<y} ■ 



(5.6) 



exp 



-A 



0- (a;. A) 



The task is therefore to find Z such that 



(5.7) 



'/'-(y,A) 



h(y) 
h'(y) 



^f(y,A)=0_(y,A)--^</)'_(y,A) 



h'{y) 



and a distribution u for Zq such that 



(5.8) 



To find the former, we start with Biane's observation that Equation (|5.7p implies 

Id 



s'iy) Ay 



So the "Riccati version" of Equation (|5.7p is 



U^{y,\) 



m{y)s{y)- 



'^wc^-{y.^) 



0_(y, A) - ^<^'_(y,A) </,_(2/, A) - ^0'_(2/, A) 

s(y) 2/a(2/)A^_(y,A) _ 2A/a(2/) 



s'(2/) </._(y. A) - i^0'_(2/, A) C/-(y, A) - ^ ' 

Table [T] then shows that, by taking Z — S'h {X)i we obtain a diffusion such that 
this does indeed hold, and hence also Equation (j5.5p — given the right v. 



6. Some deterministic examples 

There are two methods for calculating the continued fraction coefficients: the first 
uses the algorithm of [jU the second proceeds by tracking the diffusions produced 
by the recurrence 

X ^ (X) ^ o srn„{x)^ ... . 

The latter method provides greater insight but is not always applicable; the former 
method can always be used, but the determination of the constants of integration 
can be tedious. In this section, we study some well-known diffusions in deterministic 
environments which will serve to illustrate the two approaches. The processes that 
we consider have the remarkable property of belonging to a parametrised class that 
is closed under both the (/()±(-, 0)-transform and Krein duality; see Table [21 Our 
main reference is [9], Appendix I. 
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6.1. Brownian motion with drift. Let be Brownian motion with drift /n and 
set X = B^. The generator is 



1 d^ _^ ^ d 



2 dx'^ dx 

and the state space is M; the endpoints are natural and so there are no boundary 
conditions. We deduce that the Riccati equation has a unique positive solution and 
a unique negative solution. We find 



exp 



(^-fj. =F V^i^ + 2A^ and so U± = -/i =F V^i^ + 2A . 



The Stieltjes-Perron inversion formula yields 



<T±{x, dz) — dz , z > /i /2 

TT 2z 

and 

a±{x,m^\M±A ■ 

Condition ((M)) is therefore satisfied unless /i = (standard Brownian motion). In 
fact, the expansion of U± is easily obtained by elementary means: assuming for 
definiteness that /i > 0, we have 

2A 

U- = V/i^ + 2A - /i = 



2A 2A 



2/i + U- 2A 
2^+ 



2/i 



2A 



2/^- 



and 

-U+ = V/i^ + 2A + /i = 2^ + [/_ = 2/i 



2A 



2/i 



2A 



2^ 



2A 



2^ + • • • 

Let us demonstrate how these expansions can be obtained by the algorithm of 
Again, for definiteness, suppose that ^ > 0. We need to distinguish two cases, 
namely wq = and mq 7^ 0. 

Take mq = 0. Then hi = — /i, Wi — —fix and so 

where c is a constant of integration. There is no possibility of choosing this constant 
so that ui is negative. But by taking c > 0, we obtain ui > 0. Furthermore, 

2^2; /,,2 -ie^'"^^ + - + c^e-^A"^ 

uJx) = — — = ^ ^ 

/ e^M^/w^ da; gi,e2M^ + ^ ^'^"^^'^ + ^ 

where C is another constant of integration. The only way to ensure that U2 is also 
positive is by taking c — Q. Hence 



Ui{x) — 2/i and U2{x) = — 



1 ^2ux 
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By iterating, we obtain the expansion of 

The other possibihty is to take uq ^ 0. Since Wo = fJ-x, we find 

p-2fj.x 

uq{x) = 

where c is a constant of integration. No choice of c can make uq positive but, by 
taking c < 0, we can make it negative. Then 

The numerator is positive; the denominator is 

5^e-2^^ - -X + — e^'^^ + C , 

where C is another constant of integration. This expression cannot be negative for 
every x unless c — 0. We deduce 

_1 p — 2/ia; 
£-2^2; + (7 ' 



uq{x) — —2/i and wi(a;) — _^ 

By iterating, we obtain the expansion of [/_(-. 

Next, we demonstrate how the expansion may be interpreted in terms of a se- 
quence of diffusions obtained via /i-transforms and Krein duahty. We have 

0±(x,O) = exp[(-/i=F lA*l)a;] • 

Suppose that n > 0. Then (f>-{-,0) = 1, (/>+(x,0) = e-^/^^" and, from Table [H it 
is easily deduced that the (/)+(•, 0)-transform of is B'f^. On the other hand, if 
fi < 0, then (/)+(•, 0) = 1, (j)-{x,0) = e"^''^, and so it is the (/)_(•, 0)-transform of 
B'^ that yields _B~'^. Also, from Table [1] it is immediate that the Krein dual of i?^ 
is B~'^. Putting these results together, we find, for example, that the expansion of 
J7_|- for /i > corresponds to the sequence 

X = B'^ ''^+(''°)"'^""^^°''"'^ J^-fJ. Krein duality^ ^ ^ 

0,(.,O)-transfo..m^ ^_ ^ Krein duality^ ^ ^^^^ ^ 

6.2. A Bessel process. Let BES(p) denote the Bessel process with parameter p 
and, for p > 0, let X = BES(p). Then 

^ 1 d^ p + i d 
2 ax^ x ax 
is entrance-not-exit, oo is natural, 

(j)-{x, A) = x-PIp {y2\x^ , (l)+{x, A) = a;-PA'p (%/2Ax) 

and 

A) = V2A^^±ii^, -[/+(..,A) = 72A^^+^(^^). 

To obtain these expressions for the Riccati solutions, we have made use of the 
three-term recurrence relations satisfied by the Bessel functions Ip and Kp: 

2p 

Ip-i{z) - Ip+i{z) = —Ip{z), 2I' {z) = /p_i(z) + /p+i(z) 
z 



DIFFUSION PROCESSES AND CONTINUED FRACTIONS 21 



X 


Parameters 


X- 


X+ 


X* 




< 






B-f^ 




> 




B-i^' 


B-f^ 


BES(p) 


P < -1 


BES(-p) 


BES(p) 


BES(-p- 1) 


BES(p) 


-Kp < 0, 
killing 


BES(-p) 


BES(p), 
killing 


BES(-p-l), 
reflecting 


BES(p) 


-Kp < 0, 
reflecting 


BES(p), 
reflecting 


BES(p), 
reflecting 


BES(-p-l), 
killing 


BES(p) 


p = 


BES(O) 


BES(O) 


BES(-l) 


BES(p) 


<p < 1 


BES(p) 


BES(-p), 
killing 


BES(-p- 1) 


BES(p) 


p > 1 


BES(p) 


BES(-p) 


BES(-p- 1) 



Table 2. The 0±(-, 0)-transform X"^ and the Krein dual X* for 
some processes X in a deterministic environment: B^^ denotes 
Brownian motion with drift /i and BES(p) denotes a Bessel process 
of parameter p. 



and 

Kp^iiz) - Kp+iiz) = -—Kp{z) , -2K'{z) - Kp^i{z) + Kp+^{z) . 

z 

The Stieltjes measures are pT| : 

oo ^ 

cr_(a;,{0}) = 0, a_(x,dz) = ^-(5,2^ , 

fe=i '2i^ 

where the jp^k are the positive zeroes of the Bessel function Jp, and 

cr+(a;,{0}) = - , (j+{x,Az) = —^ — 



x' ^'zx J^{^x)+Y:^{V2zx) 

Since the jp^k increase linearly with fc as fc — > oo, (T-{x) satisfies the moment 
condition. On the other hand, since 

(T+ (x, dz) ^ c(p, a;) z^'^^ dz as z — , 

the moment condition is not satisfied by cr+(a;, •). By using the second of the 
recurrence relations satisfied by the Bessel functions, it is straightforward to verify 
that 

U^/2Xx) 2(p+l) /p+2(V2la:) 

V2\x Ip+i(\/2\x) 

2A 2A 



2(p±i)^^ W(7g.) ^(^^^ 2A 

^ Ip+l{\/2\x) 1 



2(p+2) 

X ^ 2(p+3) , 
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Also, despite the fact that cr+{x) does not satisfy the moment condition, 



2A- 



2p 

X 



2A 



2(P-1) 



2A 



2(p-2) 



2A 



2(p-3) 



This expansion does not contradict Theorem [1] for the coefficients cannot aU be 
positive. 

It is also possible to obtain both expansions from the algorithm of SjH First, 
consider the choice uo = 0. Then h\ = —h and 



^~2W^{x) ^ ^2p+l 



We obtain 



ui{x) 



X 



2p+l 



2.2p + 2 

2p+2 



for some constant of integration c. Then 



U2{x) 



-2Wi /„2 



lu{ 



,2p+l 



2p+2 



CX 



-2p-l 



rp-2Wl/,,2j^ 3;2p + 4 , c2 „ 

•I ^ /"lUJ- 4(p+l)2(2p+4) + 2ff2 ~ 2^ + 



where C is some other constant of integration. U2 cannot be positive unless c — 0. 
Hence 

2p+2 



ui{x) 



and 



U2{x) 



2p+i ( X y 



4(p+l)2(2p+4) 



c 



etc. 



Turning now to the case Uq 0, we have 

a;-2p-i 

uo{x) = - 



^ 2p 

for some constant of integration. Without going into details, let us simply say that 
Ml cannot be of one sign unless c = 0. By iterating this argument, we deduce the 
expansion of U+. 

Next, we examine the sequence of diffusions associated with these expansions. 
With the help of Table [2l we see that the diffusions associated with U- are: 

X = BES (p) BES ip) ^'•""'^"""^^> .%„ (X) = BES {-p - 1) 

0_(.,O)-t..ansfon.^ BES {p + 1) '^"^'"'"> O (X ) = BES {-p 2) CtC. 
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When applied to C/+, using uo,+ — —2p/x, the first iteration of the algorithm 
yields 

X = BES ip) ^+(--")~^°"^ BES i-p) {X) ^ BES (p - 1) . 

This is precisely the Ciesielski-Taylor pairing generalised by Biane [6]. 

Other well-studied diffusions that lead to simple continued fraction expansions 
are exponential Brownian motions and squared Bessel processes. 

7. Diffusion in a Brownian environment with positive drift 
Now, let a = 1 and suppose that 

W{x) = ^iX + Bx , M > : 

where B is standard Brownian motion. For every realisation of W, the process with 
generator (|1.3p is a linear diffusion and the algorithm of §4] produces a continued 
fraction whose coefficients m„(x) are random variables. We shall be interested in 
the stationary distributions of U and the u„. 

In order to study these stationary distributions, we require a stationary Fokker- 
Planck (Forward Kolmogorov) equation for the "master" system of stochastic equa- 
tions 

(7.1) ^=a(y)+i?'by, 

where y : M — > M'', a : K'' — > R*^ is a vector- valued function with components a^, 
and b is a (fixed) diagonal d x d matrix given by 

b diag(6i) . 

The flow associated with Equation (|7.ip . interpreted in the sense of Stratonovich 
[35], defines a Feller semigroup whose infinitesimal generator is given by 

Its adjoint is the Fokker-Planck operator 



(7.3) 



If / is a normalisable solution of the Fokker-Planck equation 

(7.4) ^^f = 0, 

then it is the density of a stationary solution y of the stochastic equation (|7.ip . 

7.1. The stationary distribution of the continued fraction. The equation 
satisfied by the Riccati variable is 

(7.5) ^ = 2X-2^iU -U'^ -2B'U . 

ax 

This is of the form flA} with d=l, 

a = 2A-2/xy~y^ 6 = -2 . 
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Denote by fu the density of the distribution of U. The Fokker-Planck equation is 
then 

A /_ (2A - 2^ly - y') fu + 2y^ {yfu)} = . 



dy 



dy 



This has one normahsable solution, namely 
(7.6) /c/(2/)=cy-^"iexp[-^-- 

. ^ y 

This is the density of the inverse Gaussian distribution 



y > 0. 



7.2. Proof of Theorem [3l The equation for the continued fraction coefficient uq 
is 

duQ 2 / 

— — = — 2/iUo ~ Ur. ~ 2B uq ■ 
dx 

This is the homogeneous version of Equation (j7.5p : for /i > 0, it has no normahsable 
solution. Hence uq = 0. The equation for ui is then 



(7.7) 



dui T 

= -ul + 2B'ux . 

dx 



This is of the form ([7T|) with d = 1, 

a = 2fiy ~y^ , b = 2. 

Denote by / the density of the distribution of ui. The Fokker-Planck equation is 
then 



(7.; 



_d 

dy 



(^2fiy + y^)f + 2y—{yf)\^0 



This has one normahsable solution, namely 



(7.9) 



f{y) = cy 



dy 



, y>0. 



This is the density of the gamma distribution. 

More generally, for c? e N, the stochastic equation satisfied by the first d of the 
non-zero it„ is of the form (17.11) with 



A. + £(-l)V- 



y1 and h = 2(-l)' 



Equation (|7.ip for the density, say fd^ of the joint stationary distribution of the 
random variables ui, . . . , Ud becomes 



d 



d 



We are now in a position to prove Theorem |3l Let / be defined as in Equation 
7.9p . We will show that the probability density function 



fd{yi, ...,yd) ■■= YifiVi) 
solves the Fokker-Planck equation (|7.10p . 
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Proof. We proceed by induction on d. The case d — 1 has aheady been dealt with. 
We make the induction hypothesis: namely, we set 



fdivi, ■■■,yd) ■=Wf{yi) 

i=l 

and suppose that Equation (|7.10p holds for some d. For rf + 1 , the left-hand side of 
the Fokker-Planck equation may be written as the sum 

A + B + C 

where 

(7-11) A:= ^ A / +2(-l)-i„^^(-l))-iA|„J,^,]l 



and 

d 

(7.13) C 



dvd+i y ^ % J 



First, we note that 

fd+i{yi, • ■ • , yd+i) = fdiyi, • ■ • , yd)f{yd+i) ■ 

So 

A = f{yd+i) E |- 1 -«^/^ + 2(-i)^-^?/. E(-iF^^^ [y^M I 

= f{yd+i) -0 = 

by the induction hypothesis. Next, we remark that 

B = J2 TT- {2(-i)^-'y.(-i)' [yd+ifd+i]} 

ayd+i j-[ dyi 



__d 

dy. 



— I {-l)Sd+if{yd+i) E(-l)'^'2 A [yj,] I 
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Hence 

A + B + C = 

-ad- 



^{(-i)'^..../w)|:(-ir-2AM 

d „ 

-ifiyd+i)fd + i-lfyd+ifivd+i) E(-l)'"'2 a- [y^fd] 



+2{-i)%+i{-^)''fd^^[yd+if{yd+i)] 
oyd+i 



d 



dyd^ 



d 

-ad+if{yd+i)fd + 2yd+ifdT, [yd+if{yd+i)] 

Oyd+i 



Now, 



dyt 



[yifd] 



dy. 



[yiliyi)] 



2yj^ [y*/(yj)] 



?;,(2pt - yi)f{yi) ^ -{2fi- yi)fd, 



where we have used the fact that f{yi) solves Equation (|7.8p with 1 replaced by i. 
Therefore 



A + B + C = /d 



-ad+i/(2/d+i) + 2j/d+iTi [yd+i/(«/d+i)] 

By using the definition of Ud+i, it is easy to verify that 

d 

2{-lfyd+i ^(-1)^-1(2^^ - y,) - Ud+i = yd+i{yd+i ~ 2/i) . 

Hence 

A + B + C 



fd 



d 



dyd 



+1 



2yd+] 



d 



dyd+] 



[yd+if{yd+i)] - yd+i(2M - yd+i)f{yd+i) 



/d • = 



and the proof is complete. 



□ 



If we assume that this Fokker-Planck equation admits no more than one smooth 
solution, then it follows that the m„ are, in the stationary regime, independent 
and have the same gamma distribution. The formula (|7.6p reproduces the result 
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obtained by Letac & Seshadri Since all the coefficients are positive, this 

continued fraction yields twice the Laplace exponent ip^. In particular, we can 
assert that, in a Brownian environment with positive drift, in the stationary regime, 
the reciprocal of the mean duration of the excursions below the starting point is 
gamma-distributed. The corresponding random measure cr_ was studied by Marklof 
et al. in [32] : they found that the essential spectrum is [0, cx)), with an empty 
absolutely continuous part. 
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